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Abstract 



Given the density matrix p of a bipartite quantum state, the quantum separabihty problem 
asks whether p is entangled or separable. In 2003, Gurvits showed that this problem is NP-hard 
if p is located within an inverse exponential (with respect to dimension) distance from the border 
of the set of separable quantum states. In this paper, we extend this NP-hardness to an inverse 
polynomial distance from the separable set. The result follows from a simple combination of 
works by Gurvits, loannou, and Liu. We apply our result to show (1) an immediate lower bound 
on the maximum distance between a bound entangled state and the separable set (assuming 
P 7^ NP), and (2) NP-hardness for the problem of determining whether a completely positive 
trace-preserving linear map is entanglement-breaking. 

1 Introduction 

Once dubbed "the characteristic trait of quantum mechanics" p], the phenomenon of quantum 
entanglement has been (theoretically) harnessed in a multitude of quantum computational tasks 
and areas, including quantum teleportation j2j, superdense coding [3_, quantum parallelism (e.g. 
as in Shor's factoring algorithm |1]), quantum communication complexity [5l [6], and quantum 
cryptography [7]. In response, there have been a number of proposed entanglement detection 
criteria and measures, such as the positive partial transpose (PPT) criterion [S] IS], von Neumann 
entropy of reduced states [ID], entanglement of formation [llj, relative entropy of entanglement |il2j . 
entanglement cost [11] , distillable entanglement [11] , negativity [13] , and the first need for positive 
but not completely positive maps in physics 0, to name but a few (see [TJl [151 US] surveys). 
Yet, the problem of determining whether an arbitrary quantum state is entangled or not (where in 
the latter case, the state is called separable)^ dubbed the Quantum Separability problem (QUSEP), 
was proven NP-hard by Gurvits in 2003 p7]. 

Let us discuss how one formulates QUSEP in a slightly more formal manner. Let I?(C^ C^) 
denote the set of bipartite density operators acting on C*^ , where M and N denote the 

dimensions of the respective subsystems. A quantum state p G !D(C^^ C'^) is called separable if 
and only if it can be written 



k 

for unit vectors la^) G and ^ C^, and real vector p such that Pfc > for all k and Y^i^Pk = 1- 
The latter constraint implies that the set of separable density operators, which we denote by Sm,n, 
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is a convex set, being the convex hull of all pure product states |a) |6) G (8) C'^. Intuitively, 
QUSEP is thus the problem of determining whether a given state p € P(C^ C'^) is in Sm,n- 
However, a problem arises when we encode our state p using a computer — due to the constraint 
of finite precision, we cannot in general encode the density matrix of p exactly. In particular, if p 
sits on the border of Sm,n, it may be that the slightly perturbed density matrix we actually encode 
now sits slightly outside of Sm,n- This makes the problem ill-defined. 

To circumvent this problem, one solution is to allow a margin of error in the vicinity of the 
border of Sm,n- This formulation is known as Weak Membership. Roughly, the Weak Membership 
problem over a convex set K C IR™ (denoted WMEM^(X), and defined formally in Section[2]) asks 
to decide whether a given point y G R"^ is in K, with the proviso that an algorithm is allowed to err 
on points lying within some fixed Euclidean distance /? > from the border of K. Thus, to make 
QUSEP well-defined, we consider the formulation WMEMi3{Sm,n)- Observe that Sm,n 2 IK,™ 
required for WMEM^(i^) — we deal with this explicitly in Section |2] by more correctly letting 
Sm,n denote the set of real Bloch vectors [18j corresponding to the elements of I?(C*^ (gi C^). 

In 2003, Gurvits showed [17^ that WMEM/3(5M,Af) is NP-hard via a polynomial time Turing 
reduction from the NP-complete problem PARTITION. Intuitively, a Turing reduction describes 
how to solve a problem A (e.g. PARTITION) by running an algorithm for a second problem B 
(e.g. WMEM/3(5M,Af)) possibly multiple times. PARTITION is defined as the problem of deciding 
whether a finite set of integers can be partitioned into two sets of equal sum. However, PARTITION 
is known to be NP-hard only if the magnitudes of the input integers are exponentially large with 
respect to input length — otherwise, the problem can be solved efficiently using a dynamic program- 
ming approach [W. It follows, as observed by Aaronson and later documented by loannou [20j, that 
the reduction of Ref. [17J shows NP-hardness for WMEM/3(5Af,Ar) only when /? < l/exp(M,iV), 
i.e. when the input state is allowed to be exponentially close to the border of Sm,n- 

In an attempt to strengthen this result, Gurvits then devised (as explained in [20]; see also j21j ) 
the following reduction from the NP-complete problem CLIQUE (defined in Section |2]): 

CLIQUE <m RSDF WVAL„(5Af,7v) <t WMEM^{Sm,n)- (2) 

Here, RSDF is the problem Robust Semidefinite Feasibility (defined in Section [2|, 'WYALa{SM,N) 
is the problem Weak Validity [22] (which intuitively asks one to decide whether a given hyperplane 
is a separating hyperplane for a given convex set modulo some error a > 0, and whose precise 
definition is not needed here), <t denotes a Turing reduction, and <m denotes a many-one re- 
duction. A many-one reduction is a special case of a Turing reduction in which the algorithm 
for problem B is invoked only once, the output of which is immediately returned as the output 
for problem A. Unfortunately, the link WVALa{SM,N) I^T WMEM/3(5A.f,Af) above is based on the 
Yudin-Nemirovskii theorem f23], which uses the shallow-cut ellipsoid method, and also results in 
exponential scaling for f3 ||20l. Thus, this reduction again shows NP-hardness of WMEM^(5j\/,Ar) 
only for P < l/exp(M,7V). 

The main result we show in this article is as follows. 

Theorem 1. WMEM/3{Sm,n) is NP-hard for (3 < l/poly(Af, A^), i.e. is strongly NP-hard. 

A problem is called strongly NP-hard if it is NP-hard even if the magnitudes of its numerical 
parameters are polynomially bounded in the length of its input [H] . To show Theorem [l| our obser- 
vation is that we can replace the component 'WYALa{SM,N) <T WMEM/3(5m,7v) in Eq. Q with a 
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recent non-ellipsoidal Turing reduction of Liu [21] (see also Ref. [25], as discussed in Ref. |24j^) from 
the problem Weak Optimization (denoted WOPTe(i^), and defined in Section [2]) to WMEMp^K). 
We thus use the new reduction chain: 

CLIQUE <m RSDF <m WOPT,(5m,7v) <t WMEM^(5m,7v). (3) 

To make this work, our technical contribution is the reduction RSDF <m WOPT^{Sm,n), which 
uses ideas similar to those in the reduction RSDF <m WVALa{SM,N) [20] . 

This article is organized as follows. In Section |2] we introduce all necessary definitions and 
notation. Section [3] presents the proof of Theorem [T] In Section [4] we discuss two applications 
of Theorem [T| We first apply the positive partial transpose (PPT) entanglement detection crite- 
rion [H [9] to obtain immediate lower bounds on the maximum Euclidean distance between a bound 
entangled [2F state and Sm,n (assuming P 7^ NP). We next use the Jamiolkowski isomorphism [27] 
to show NP-hardness of the problem of determining whether a completely positive trace-preserving 
linear map (i.e. a quantum channel) is entanglement-breaking [28]. We conclude in Section [5j 

2 Definitions and Notation 

In this section, we formally define the following problems needed to show Theorem [TJ CLIQUE, 
Robust Semidefinite Feasibility (RSDF), Weak Optimization (WOPTe(-fC)), and Weak Membership 
(WMEM^(ir)). All norms are taken as the Euclidean norm || II2 (where ||^||2 corresponds to the 
Frobenius norm if ^ is a matrix). The letter Q indicates the rational numbers. The notation := 
is used to indicate a definition. We denote (column) vector v by v, its conjugate transpose as v''', 
and its zth entry as Vi. We use the notation {a) to signify the number of bits necessary to encode 
an entity a. Specifically, if a = a/b is rational, we define (a) = (a) -|- (6), and for matrix A, we let 
(A) = "^-j (Aij) (similarly for vectors). 

First, the NP-complete problem CLIQUE is stated as follows. 

Definition 1 (CLIQUE). Given a simple graph G on n vertices, and c < n, for n,c £ decide, 
with respect to the complexity measure (G) -|- (c): 

If the number of vertices in the largest complete subgraph of G is at least c, output YES. 
Otherwise, output NO. 

Here, we take (G) = (Ac), where Aq is the n x n adjacency matrix for G, such that ^^[2, j] = 1 if 
vertices i and j are connected by an edge, and ^^[2,7] = otherwise. Next, we have the problem 
Robust Semidefinite Feasibility. 

Definition 2 (Robust Semidefinite Feasibility (RSDF)). Given k rational, symmetric / x / matrices 
Bi,...,Bk, and ^ Q, with C,r? > 0, define g{Bi,...,Bk) := maXxg]ai,||x||2=i E!Li(x'^^ix)^. 
Then, decide, with respect to the complexity measure Ik + Yli=i (Bi) + (C) + iv)'- 

If g{Bi, ...,Bk)>C + V, output YES. 
If g{Bi, ...,Bk)<C-V, output NO. 

^Both Ref. [24] and Ref. [25] only pertain to the third reduction in Eq. ([3| — i.e. they discuss the reduction 
WOPT,(A') <T WMEM^(7^) for arbitrary K. 
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Figure 1: An intuitive picture of the sets S{K,5) and S{K,—6), respectively. We use the term 
"—5" in the latter illustration to stress the fact that K is contracted. 

We have assumed C ^ without loss of generality above, since g{Bi, . . . , Bk) > 0. This will be 
necessary later in Lemma |4] when we need to take g{Bi, . . . , B^). We have also defined RSDF as 
a promise problem, meaning we are promised the input will fall into one of two disjoint cases which 
may be separated by a non-zero gap, and we must distinguish between the two cases. One could 
equivalently lift the promise and allow input falling in the "gap" region (i.e. C~V < 9{Bi, ■ ■ ■ , -Bfc) < 
C + v) — foi" such input, we would consider any output to be correct (i.e. YES or NO). 

Moving on, in order to discuss the Weak Optimization and Weak Membership problems, we 
first require the following definitions. Let K C R*" be a convex and compact set, and define the 
sets S{K, (5) := {x G I 3 y E s.t. ||x - y < 5}, and S{K, -6) := {x e K | S(x, 6) C K}. 
Roughly, S{K, 5) can be thought of as extending the border of X by 5 in Euclidean distance, and 
S{K, —6) can be thought of as taking the core of K, which is 5 away from the border of K. The two 
sets are depicted in Figure [T} As per Ref. p2], we shall require that be a well-bounded p-centered 
set, meaning that K C S{0, R) for G the origin and some rational R > 0, and ^(p, r) K for 
a known point p G X and some rational r > 0. This ensures K is full-dimensional and bounded. 

Finally, we set the encoding size of K as {K) = m + {r) + (R) + (p) jl7j . The presence of 
parameter m in this expression allows us to perform vector addition and scalar multiplication on 
elements of K in polynomial time. We can now formally define Weak Optimization and Weak 
Membership over the set K as follows. 

Definition 3 (Weak Optimization (WOPTe(i^))^). Let K C R™- be a convex, compact, and 
well-bounded p-centered set. Then, given c G with ||c||2 = 1, and 7,e G Q, such that e > 0, 
decide, with respect to the complexity measure (K) + (c) -|- (7) -|- (e): 

If there exists y G S{K, — e) with c^y > 7 -|- e, then output YES. 
If for all x G S{K, e), c'^x < 7 - e, then output NO. 

Definition 4 (Weak Membership (WMEM^(K))). Let K C R"* be a convex, compact, and well- 
bounded p-centered set. Then, given y G Q"*, and error parameter /? G Q, such that /? > 0, decide, 
with respect to the complexity measure (K) + (y) -|- (/?): 

If y G S{K, -/?), then output YES. 
If y sIk, /?), then output NO. 

Both of these problems are also stated as promise problems. Roughly speaking, the former asks 
whether there exists a point in the "core" of K that achieves a threshold value slightly higher than 
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7 for a linear function defined by c. The latter asks to distinguish whether a given point y is in 
the "core" of K or "far away" from K. We remark that unlike here, in Ref. [21] the inputs to 
WOPTe(i(r) and WMEM^(i(') are real (as opposed to rational), and specified using poly(m) bits of 
precision. The latter is because such precision suffices if one demands e and (5 to be at least inverse 
polynomial in the input size [23], i.e. if one allows at least "moderate" error (which we will also 
demand here). It is easy to see that we can exactly represent any such poly(m)-bit real numbers 
as rational numbers in poly-time using poly(m) bits as well, and hence the case of Ref. |24j is a 
special case of our definitions here. 

With the definitions above in place, our goal in Sectionjsjis to show the reduction of Eq. (|3]), i.e. 
that an instance of CLIQUE can be solved in polynomial time given an algorithm for WMEM^(5A/,Ar). 
Before proceeding, it remains for us to verify that the set of separable quantum states in D(C*^(E>C^) 
satisfies our requirements on if C K,™ in the definitions of WOPTe(-fC) and WMEM^(if). To do 
so, we first represent our quantum states as real vectors in R"*, for m set as follows. Let Hm,n 
denote the set of Hermitian operators mapping C*^ (8) C''^ i— > C'^ (Xi C^. Then, via the isomorphism 
between Hm,n and R^'^', we can write any p £ V{C^ C^) as (e.g. [T8]): 

1=1 

where we have chosen as a basis for Hm,n the identity and the traceless Hermitian generators of 
SU{MN), the latter denoted by ai and such that Tr((Ti(Tj) = 26ij (for 6ij the Kronecker delta) [T8] . 
The vector r G IR,^^ ^ ~^ is called the Block vector of p, whose terms are given by rj = Tr(po"i). 
Henceforth when referring to Sm,n^ we shall mean the convex set of Bloch vectors in E,^ ^ ~^ 
corresponding to separable states in V{(C^ (8) C-^). Having represented Sm,n in terms of real 
vectors, we require it to be p-centered and well-bounded. This follows since Sm.n is contained in 
an origin-centered ball of radius R = y^2{MN — 1)/MN [16^ , and contains an origin-centered ball 
of radius r = ^J2/MN{MN - 1) p9] , where the origin E R*^ ^ ^ corresponds to the maximally 
mixed state, and R and r are with respect to the Euclidean norm. We remark that the extra factor 
of \/2 appearing in the expressions for r and R is due to switching from the density matrix to 
the Bloch vector representation (see Lemma [9] in Appendix |A| — this scaling will not affect our 
analysis. We further require that Sm,n be compact, which holds since the set of pure product 
states is closed and bounded, and the convex hull of a convex compact set is also compact [30]. 

Finally, for concreteness, we remark that from the values of m, r, R, and p above for Sm,n, it 
follows that for the definitions of W0PT^(5M,Af) and WMEMf3{SM,N) (Def . [s] and |4) respectively), 
{Sm,n) = m + {R) + (r) + (p) < poly(MA^). ' 



3 The Reduction 

We now show Theorem [T] by demonstrating the polynomial-time Turing reduction of Eq. ^ . Since 
CLIQUE is NP-complete, this will imply NP-hardness of WMEM^(5A/,Ar). In addition, we must 
show that in Eq. ^ , one can choose e, /? > 1/ poly(M, N) — the former is required for the reduction 
to run in polynomial time (due to the run-time of Theorem [5|, and the latter is required to prove 
strong NP-hardness (since the case of /3 < l/exp(M,N) is already known to be NP-hard [E]). We 
proceed by considering each link of Eq. ([s]) in order. As will be discussed, the ffi'st and last links 
are provided by Ref. [20] and Ref. |24j , respectively. Our job is to show the second link. 



5 



The first link in Eq. (|3j) is given to us by the following theorem |20j . Unless otherwise stated, 
by a poly-time reduction, we mean with respect to the encoding size of the problem instance, as 
defined in Section [2| We use the notation 11 = {input parameters) to denote an instance 11 of a 
given problem, with 11 specified by the given input parameters. 

Theorem 2 (Gurvits and loannou [IQ]). There exists a poly-time many-one reduction which maps 
instance Hi = {G,n,c) o/ CLIQUE to instance II2 = {k,l, Bi, . . . , B^XtV) 0/ RSDF, such that 
k = n{n- l)/2, l = n, Bi^ Q"''" and \\Bi\\^ G 9(1) for aUl<i<k, ( = 9(1), r] G n{n^'^). 

We refer the reader to Ref. [2D] for the details of the proof of this theorem, but highlight that it 
relies heavily on the following theorem of Motzkin and Straus ||3l] , which ties the maximum clique 
number of a graph to optimization of a square-free quadratic form over the unit simplex: 

Theorem 3 (Motzkin and Straus [31j). Denote by {i,j) £ G an edge in graph G between vertices 
i and j, and let uj be the order of the maximal complete graph contained in G. Let A„ denote the 
simplex A„ := {x G E." | > 0, ||x||j^ = 1}. Then 

max XiXj = - I 1 I . (5) 

For later reference, we remark that the matrices Bi G R"^" from Theorem [2] will have the following 
structure — to each Bi, we uniquely assign an index (s,t) from the adjacency matrix Aq of G, 
1 < s < i < n, such that Bi has all entries zero, except for entries (s, t) and (t, s), which are set to 
the (s,t)th entry of Aq- We hence require k = n(n — l)/2 matrices -Bj, as claimed by Theorem |2j 
We now demonstrate the second link in Eq. ([3|. 

Lemma 4. There exists a poly-time many-one reduction which maps instance Hi = {k,l, Bi, . . . , B^, C) 
of RSDF to instance II2 = (c, 7) of WOPT^{Sm,n), where we define for convenience A := 

2X^^=1 Il-^ill2' '^^^ such that: 

> M = k + l 

• c = c/ ||c||2 for some c G Q"^ with ||c||2 G 0(m-^/^A) and m = M'^N'^ — 1 

' 7 = ^(VCT^ + \/C^) 



4||c|l2{AW-l)+l 

Proof. The heart of the mapping from Hi to II2 is given in Refs. |17| |20] . and involves rephrasing 
the function g{Bi, . . . ,Bk) from RSDF in terms of convex optimization over the set of separable 
density matrices acting on C*^ (8) C^. We briefly summarize this here for later reference. Let 
M = k+1, N = 1, and consider the matrix G G E*^^""*^^, such that 



G :-- 



( Ai ... Am-1 \ 
Ai ... 

V Am-1 . . . ) 



(6) 
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where each Ai E R^^-^ is symmetric and all zeroes except for its upper-left I x /-dimensional 
submatrix, which we set to Bi. It is easy to see that ||C||2 = A, as defined in the statement of our 
claim. One can then write (Proposition 6.5 in Ref. [T7]): 



■\/g{Bi, ...,Bk) = max 



M-l 



(xT^,x)2 = rnax TV(C7/,sep), (7) 



xeR^,||x||2=i \ ^ pscpefCcA-fcac^) 

where psep denotes a separable density matrix. Thus, Hi is reduced to optimizing the linear 
objective function Tr(C/o) over all separable density matrices psep G P(C*^ C^). This concludes 
the referenced work of |17| |20] . 

What remains is to explicitly rephrase the problem in terms of Bloch vectors and apply simple 
convex geometric arguments to complete the reduction, as well as characterize scaling of the error 
parameter e. To do so, first use Eq. Q and the fact that Tr(C) = to write: 

/ / a-Pn^-i \ \ M^N^-l 

Tr(Cp)=Tr (c( — + - ''■'^01=2 ^ H • Tr(Ca,) = e^r, (8) 

for (jj the generators of SU{MN), Ci := ^Tr(C(Tj), and r the Bloch vector of p. Set m = M'^N^ — 1, 
c = c/ llcllg. In terms of Bloch vectors, our objective function Tr(Cp) in Eq. ([T]) can hence be 
rephrased as /(r) := c'^r, with /max := ^SiXreSm n /('^)- We remark that unless C is the zero matrix 
(i.e. the input graph to CLIQUE has no edges), we have ||c||2 > 0. Also, since Tr((TjCJj) = 26ij, it 
follows from Eq. (|8| and the Cauchy-Schwarz inequality that ||c||2 G ©(m^/^A). 

To complete the reduction, we now must show the following (for 7 and e to be chosen as needed): 
If /max > l|c||2 ^ VC + Vi then there exists an r G S{Sm,n, — e) such that /(r) > 7 -|- e (i.e. a YES 
instance of RSDF implies a YES instance of WOPT^{Sm,n))- If /max < l|c||2 ^ VC ~ then for all 
r G S{Sm,n, e)) /(r) < 7 — e (i-e. a NO instance of RSDF implies a NO instance of WOPT^{Sm,n))- 
The ||c||2 ^ term in these expressions follows from our definition of c, and the square root in 
follows from the square root in Eq. (7). We proceed case by case. Set 7 = 2||e|| iVC + V + VC ~ v)} 
and let us choose e as needed. 

• Case 1: /max > pjj^VC + V- 

Let r* G Sm,n be such that /(r*) = /max- To find an r G S{Sm,n, — e) such that /(r) > j + e, 
we first use the fact that for any well-bounded origin-centered convex set K C IR'", it holds 
that for all x ^ K, there exists a y G S{K, — e) such that ||x — y||2 < 2eR/r [22] (where R 
and r are the radii of the ball containing K and the origin-centered ball contained within K, 
respectively). Plugging in the definitions of r and R for Sm,n from Section |2] it follows that 
there exists an r G S{Sm,n, — e) such that ||r — r*||2 < 2{MN — l)e. Since / is linear, we can 
write: 

|/(r) - /(r*)| = |cT(r - r*)| < ||c||2 ||r - r*||2 < 2[MN - l)e, (9) 

where the first inequality follows from the Cauchy-Schwarz inequality. Thus, in order to have 
/(r) > 7 + e as desired, it suffices to have 

/(r) > /max - 2{MN - l)e > 7 + e, (10) 
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into which substitution of our values for 7 and /max gives that setting 



^ - 4||c||2(MAr- 1) + 1 ^ ' 



suffices to conclude we have a YES instance of W0PTe(5 



Case 2: /max < pjj^VC - V- 

We must show that for an appropriate choice of e, we have /(r) < 7 — e for all r E S{Sm,n, e)- 
Set r G S{Sm,n,£)- Then by the definition of S{Sm,n,£), there exists some r' G Sm.n such 
that ||r - r'llg < e. By Eq. @, it follows that 

|/(r)-/(r')| <e. (12) 



In the worst case, one has /(r') = /max- Hence, by Eq. (12), we have that /(r) < /max + e 
for any r G S{Sm,n, e)- To achieve /(r) < 7 — e then, set /(r) < /max + e < 7 — into which 
substitution of our values for 7 and /max yields that choosing 



2 ||c||2 + 2 

suffices to conclude we have a NO instance of WOFT^{Sm,n 



e < ^,7^^ (13) 



□ 



Observe that combining Theorem [2] and Lemma ^ gives M = N = "^"^"^^ + 1 • Following an 
argument of loannou (Section 2.2.5 of Ref. ^20]), one can likewise show Lemma |4] for M > N 
by padding the matrix C from its proof with extra x A^-dimensional zero matrices. Thus, the 
hardness result we show for WMEM^(5Af^Ar) by building on this link will be valid for M > N. 

The last link of Eq. (|3j) is given by applying the non-ellipsoidal Turing reduction of Liu j24] . 
which holds for an arbitrary p-centered well-bounded compact convex set K C IR™. 

Theorem 5 (Proposition 2.8 of Ref. [24). Let K C IR,™ he a convex, compact, and well-hounded 
p-centered set with associated radii {R,r). Given an instance II = {K,c,j,e) ofWOPT^(K), with 
< e < 1, there exists an algorithm which runs in time poly{{K) ,R, [1/e]), and solves H using an 
oracle for WMEM^(i^") with (3 = r^e^ /[2^'^3^m^R'^{R + r)]. 

We briefly review some of the ideas behind the proof |2l] of Theorem [5| which builds on results of 
Grotschel etal^. WOPT^(K) is first reduced to a problem called Weak Separation (WSEP,,(K)) 
over K, which roughly asks one to determine whether a point p G IR*" is approximately in K, and if 
not, to return a hyperplane approximately separating p from K. This first reduction is achieved in 
two steps: First, use the given oracle for WSEP,y(i<r) to construct a weak separation oracle for the 
set K' := £ K \ c^y > 7}. Second, use the latter oracle in an iterative approach where, in each 
iteration, we check if a candidate point p is in K' , and if not, we update p using the separating 
hyperplane returned by the oracle. WSEP i,(ivr) is then reduced to WMEM^(i^) via three further 
reductions, the latter two of which are given by Lemmas 4.3.3 and 4.3.4 of Ref. [22]. Lemma 4.3.3, 
in particular, demonstrates that a single call to a Weak Membership oracle with two-sided error 
(i.e. the kind we use in this article) can be used to solve a variant of Weak Membership with 
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only one-sided error, i.e. where in the NO case of Def. |4j one affirms that y ^ if, as opposed to 
y^S{K,i3). 

We now show that composing the three reductions above yields a polynomial time reduction 
from CLIQUE to WMEM/3(5M,Ar), such that f3 > l/poly(M, iV). First, observe the dependence 
on [1/e] in the runtime stated in Theorem [5} We thus must be able to choose e > 1/ poly(m) > 
I/poly(M, A'^) in order for our reduction chain of Eq. ^ to run in polynomial time, where m = 
M'^N'^ — 1 for Sm,n- Let us show that we can solve our instance of CLIQUE with such a choice of 
e. Specifically, by Lemma |4j we can set 



4||c||2(MiV-l) + l' 



Piecing together Lemma [2] and Lemma |4j we first immediately have ( G ©(1) and rj G ^{1/N). 
It follows that for some positive constants ci, C2, A'^i, and we have: 



VC + ^-VC^ > aM + ]^-V^i-]^ ViV>max(iVi,iV2) (15) 



(ci + f ) - (ci - f ) 



(16) 



N 



- N{v^r+^+v^y ^^^^ 



Hence, VC + V — VC — V ^ 0,{1/N). With a little thought, we also have that ||c||2 G 0{\/ N) (see 
Appendix [A] Lemma |9]). Plugging these bounds into Eq. ( |14| yields that we can solve an instance 
of CLIQUE with e G ^{M~^N~^/'^). Combining Theorem |2j Lemma [4j and Theorem [Hj we thus 
have a polynomial time Turing reduction from CLIQUE to WMEM^(5M,Ar)- 

To show that this implies strong NP-hardness of WMEM^(5A/,Af)) it remains to ensure that 
[i > 1/ poly{M, N). Substituting our values for r, R, and m for Sm,n from Section |2] into the 
expression for /3 in Theorem [s] we have /3 = poly(M~^, iV^^, e). By our choice of e above, we thus 
have /? > 1/ poly(M, A^), as required. We hence conclude: 

Theorem 6. Given instance H = {G,n,c) 0/ CLIQUE, there exists an algorithm which solves H 
in time poly(n) using an oracle for WMEM/^{Sm,n) with parameters M = N = ""^"^ + 1 
some (5 G 0(n~'^^). More generally, for N = "^"^"^^ + 1 ^'^Z/ choice of M > N , the result holds 
for some (5 G n{M-^^N-'^^-^). 

Theorem 1 {revised). WMEM^(5M,Ar) is NP-hard for (3 < ^o\y{M-^^N-'^^'>) and M > N, or 
equivalently, is strongly NP-hard. 

We stress the phrase ^^some (3 G $7(71^''^)" in the statement of Theorem [g] above — specifically, 
we cannot have /3 G ^^(1) in our reduction, due, for example, to the expression for e in Lemma |4j 
The question of NP-hardness of WMEM^(57v/,7v) for (3 G ^^(1) thus remains open. We remark that 
the major contributor to the large negative exponent on n in the estimate for (3 is the reduction of 
Theorem [S] 
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4 Applications 



We propose two applications of Theorem [T] both of which to our knowledge were not previously 
known. In this section, for simplicity of exposition we revert back to letting Sm,n denote the set 
of density matrices corresponding to separable quantum states in I?(C^ (8> C^). By Lemma [o] in 
Appendix A, distances in the respective spaces are asymptotically equivalent. Our first application, 
which benefits directly from the improved NP-hardness results of Theorem [T] is an immediate 
lower bound on the maximum distance a bound entangled state can have from Sm,n, assuming 
P 7^ NP. To see this, recall that bound entangled states are mixed entangled quantum states from 
which no pure (state) entanglement can be distilled pGl, and are the only entangled states whose 
entanglement is capable of escaping detection by the Peres-Horodecki Positive Partial Transpose 
(PPT) |8l |9] criteriorl^ Theorem [l] implies that, unless P = NP, any test of membership for 
Sm,n must be unable to efficiently resolve Sm,n within distance P G Q{M^^^N^'^^-^) of its border 
in the general case. It follows that unless P = NP, there must exist bound entangled state(s) 
Pbe e T>{(C^ (g) C^) such that for any separable state psep S Sm,n, \\Pbe - Psep\\2 ^ f^(M~^'^A^~20-5) 
— if not, one could determine the separability of any quantum state within this region efficiently 
using the PPT test, contradicting Theorem [T| Further improvements to our hardness estimate for 
P would directly benefit this application. 

Our second application is to show that the problem of determining whether a completely positive 
(CP) trace-preserving (TP) linear map (i.e. a quantum channel) is entanglement-breaking (EB) is 
NP-hard. We remark that here we do not use the improved NP-hardness results of Theorem [T] 
(i.e. NP-hardness of WMEM^(5Af^Ar) for inverse exponential /3 would suffice for our proof below). 
At the end of this section, we briefly discuss how the improved hardness bounds of Theorem [l] 
may help extend the result here to strong NP-hardness of determining whether a channel is EB. 
Returning to our discussion, quantum channels are of interest, as they correspond to physically 
realizable processes. EB channels in particular are well-studied [321 [331 EH ISSl |36l EZ], and are 
intuitively defined as the set of channels whose action on half of any bipartite state always results 
in a separable state, i.e. they break entanglement across the bipartition. Our approach is to show 
that the ability to determine whether a quantum channel is EB suffices to solve QUSEP, implying 
the former must be NP-hard. 

We begin by stating the key ingredients required to study our problem, and then sketch a proof 
of the desired result. To simplify our discussion, we do not use Weak Membership formulations 
here — the use of such formulations to make our argument rigorous is described in the ensuing 
discussion. Let L{X) denote the set of linear operators mapping a complex Euclidean space X to 
itself. Then, EB channels are defined as follows. 

Definition 5 (Horodecki, Shor, and Ruskai [2H])- A quantum channel $ : L(C*^') L(C*^) is 
entanglement-breaking (EB) if (<I> (g) I){p) G Sm,n for ah p £ L>(C*^' C^). 

Next, informally, let EBP denote the problem of determining whether a channel ^ given as input 
is EB. How do we encode ^ as input? For this, we choos^the Jamiolkowski representation ^27^ of 

^The converse question of whether all bound entangled states escape detection by the PPT test is, however, a 
major open question [H]. 

^Ahhough we represent our channels here in Jamiolkowski form, it is straightforward to move to another repre- 
sentation, such as the Kraus representation. For example, the canonical way to determine a Kraus representation 
for $ given its Jamiolkowski representation J{^) is to reshuffle the eigenvectors of J($) into matrices (i.e. Kraus 
operators) Ki. We refer the reader to [16] for further details. 
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maps, since it known that this representation provides a direct hnk to the separabihty of density 
operators [28;, as will be stated shortly in Theorem [Tj The Jamiolkowski representation is defined as 
follows. For any <I> : L(C^) i-^ L(C^^), there exists a unique operator, denoted J($) G L(C^(8)C^), 
obtained by applying the Jamiolkowski isomorphism to the action of which is given by 

J(cl>):=[$®/](|0+)(<A+|), (18) 

where |(/>"'") is the maximally entangled state \(t)'^) = Ylk=o \k)^\k) and {\k)}^~f^ is an arbitrary 

orthonormal basis for C^. If and only if $ is CP, then is positive semidefinite f38], and if and 
only if $ is TP, then Tr^(J(<I>)) = I /N [16]. Here, A and B denote the respective subsystems of 
J($). 

Let us now state the link between the Jamiolkowski representation for quantum channels and 
QUSEP, which will be the starting point for showing our result. 

Theorem 7 (Horodecki, Shor, and Ruskai [28j). Given a linear map ^ : L(C^) i-^ L(C^'^) that is 
CP and TP, i.e. J(^) is positive semidefinite and Tryi(J($)) = I /N , respectively, the following are 
equivalent: 

1. ^ is entanglement-breaking. 

2. J[^) is separable, i.e. J($) G Sm,n- 

This theorem links QUSEP and EBP in the following (immediate) way: The problem of determining 
whether a state p G !D(C*^ (8> C^) with Tia{p) = I /N is separable is equivalent to the problem of 
determining whether the corresponding channel such that J($) = /O, is EB. Thus, if the former 
problem is NP-hard, so is the latter. To show that EBP is NP-hard, we are hence reduced to 
answering the question: Is QUSEP still NP-hard if one is promised that the input state p satisfies 
Tr^(/9) = I/N? Let us remark that, in contrast, the similar problem of determining whether a CP, 
but not necessarily TP, map $ is EB is trivially strongly NP-hard by Theorems [6] and [TJ since by 
the properties of the Jamiolkowski representation discussed, dropping the TP constraint on $ in 
Theorem [T] corresponds to dropping the constraint that Tr^(p) = I/N. 

We now sketch a proof that QUSEP remains NP-hard even if one is promised the input state 
p satisfies Ti^ip) = I/N. Our approach is to demonstrate a poly-time many-one reduction to this 
problem from general QUSEP. Specifically, we show how to efficiently map an arbitrary input 
p G P(C^ C^) for QUSEP to some p' G ^(C^ C^ C^), such that TrA'(p') = I/N (where 
A' denotes subsystem ® C*^), and p is separable if and only if p' is separable across the A' /B 
split. The mapping proceeds as follows. Assume we have an efficient algorithm Q for determining 
separability in the restricted case of Tr/i(p) = I /N . Then (each step below will be explained in the 
ensuing discussion): 

1. Define the TPCP map $ : C*^ C^ O C^ O C^ such that 

$(a) := {l-p)\{)){{)\®a+p\l){l\®I/N, 

for p = 1 - 1 /N. 

2. Define the CP, but not necessarily TP, map T : ® C*^ (g) C*^ (g) (omitting 
normalization for simplicity) 

T(c7) = (/ a^^/^) a[l® ae'^'^) , 
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where fx^ = Try^/ (cr), A' and B denote the subsystems €^ and C^, respectively, and 
the identity / acts on A' . 

3. Call Q with input T($(/>)), and return Q's answer. 

To explain this reduction, we begin with Step 2. Our goal is to transform p into a state p' 
with a maximally mixed subsystem, while preserving separability. One way to directly achieve this 
is to apply the local operations and classical communication (LOCC) map T to p, which locally 
applies the inverse of Tta_{p), ensuring Tr^(T(p)) = I/N (in this case T would act analogously 
on C^). As an aside, observe that while T is CP, it is not necessarily TP. Hence, T is a 
probabilistic map (i.e. succeeds with some non-zero probability, and for this reason may be classified 
as a stochastic LOCC map, or SLOCC), and in particular is an example of a local filter, the latter 
being first considered in [39] . To see that T preserves separability, note that T can be inverted with 
non-zero probability, and so T(p) is entangled if and only if p is, since LOCC operations cannot 
create entanglement from scratch. 

However, we cannot always begin by applying Step 2, because in general Tia{p) is not full 
rank, and even if it is, it must have a low condition number in order to allow us to compute its 
inverse reliably when applying T. Here, we define the condition number |l0] of a matrix C as 
k{C) = \\C\\ ||C~^||, where ||-|| denotes the operator norm. To address this, we first preprocess p 
in Step 1 by applying which consists of two main components: mixing with the identity, and 
attaching a "marker" ancilla. The first component, mixing with the identity, ensures that not only 
is pb '■= Tryi/($(p)) full rank, but also that k{pb) ^ 3 for all N >2. To see this, observe that since 
PB is positive semidefinite, one has k{pb) = Xma.x{pB)/>^min{pB), where XmaAps) and Ainin(pB) 
denote the largest and smallest eigenvalues oi ps respectively. Combining this with the fact that 
Pb = {I — p)TrAP + pi /N yields 

^ _ (l-p)Amax(TrAp)+p/iV , (1 - p) + p/iV 

^^PB) - (1 _ p) A^i,(TrAp) + p/N - p/N ' ^'^^ 
where we have used the trivial bounds Aniax(Tryip) < 1 and Ainin(Tr^/9) > 0. For p = 1 — 1/N, 



the right side of Eq. (19) equals (2A^ — 1)/(A^ — 1), which is upper bounded by 3 for > 2, as 
desired. The second component of the "marker" ancilla, serves to preserve separability in a 
fashion similar to T by ensuring the LOCC operation $ can be inverted with non-zero probability 
(which is achieved by performing a von Neumann measurement on the ancilla in the computational 
basis). We thus have that $(p) is entangled across the A'/B split if and only if p is entangled 
across the A/B split. 

Having transformed p into the state p' := T{^{p)) which satisfies our requirements for Q, we 
finally proceed to Step 3, in which we plug p' into Q and return Q's answer. Thus, the existence 
of an efficient algorithm Q implies an efficient algorithm for QUSEP, yielding that determining 
separability of quantum states with a mixed subsystem is NP-hard. By Theorem [TJ this in turn 
implies our desired result, that EBP is NP-hard. 

A few remarks are in order. First, to make this argument rigorous, we must use Weak Mem- 
bership formulations, for which we will need to introduce additional notation. Let J^{C^ (8) C'^) := 
{p G V{C^ ® C^) I Tta{p) = I/N}, or equivalently, the convex set of operators which are the 
Jamiolkowski representations of channels $ : L(C^) i-^ L(C^'^), and let gM,N ■= (g) C^) n 

Sm,n, or equivalently, the convex subset of EB channels in .F(C*^ C^) (convexity of gM,N is 
clear from Def. |5]). In a rigorous argument, instead of reducing from QUSEP to EBP as done 
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above, one must more acrarately reduce WMEM^(5Af^Ar) to WMEM^/(^j\,/^7v) (for the same rea- 
son regarding finite precision as in Section |3]) . As per Definition |4| this means we require that 
relative to ^(C^^ C^), Qm,n is be compact, well-bounded and p-centered. To see that Gm,n 
is well-bounded and p-centered, observe that it follows directly from the fact that Sm,n is con- 
tained in a ball of radius R = a/ {MN — 1) /MN and contains a ball centered at I/MN of radius 
r = y^l/MN{MN -1) (see Section |2[ that Gm,n must also be contained in a ball of the same ra- 
dius R and contain a baig of the same radius r. To show compactness relative to ^(C*^®C^), we 
first have by the continuity and linearity of the partial trace that the complement of ^(C*^ (8> C-^) 
is open, implying ^(C*^ C^) is closed. Since Sm,n is compact and Gm,n = ^ C^) n^M.Af, 

it follows that Qm,n is compact relative to P(C^ ^^)) since in a metric space the intersection 
of a closed set and a compact set is compact (Corollary of Theorem 2.35 in [H]). Finally, since 
Qm,n Q -^(C^ C^) C V{C^^ C^), we have that gM,N is compact relative to ® C^) 

(Theorem 2.33 in [41 ). Hence, Qm,n satisfies our requirements for Weak Membership. 

To address the next technical detail, recall that Weak Membership takes as input real vectors, 
not operators. As discussed above, to ensure that Qm,n is compact, well-bounded and p-centered, 
we wish to phrase EBP relative to J-{G^ as opposed to P(C''^(8>C'^). To do so, we must map 
elements of :r(C^^ C^) to real vectors in rC^v/^-i)^', rather than E^'^'-i, where (M^ - l)N'^ 
and M^iV^ _ i the dimensions of ® C^) and I?(C*^ (g) C^), respectively (to see that the 

dimension of J^(C*^ ® C^) is (M^ - l)iV^, observe that the partial trace constraint TrA(p) = I/N 
effectively places N'^ — 1 linear constraints on /) G 2?(C*^ ® C^)). To achieve this mapping, we use 
what is sometimes called the Fano form [32] for any p G P(C^^ C^): 



MN 



1=0 j=0 



(20) 



Here, and o"^ denote the real (M^ — 1) -dimensional Bloch vector for subsystem A and vector 
of generators for SU{M), respectively, as defined in Eq. ([8|, and each real Tij can be defined 



Tij := Tr(cj^ (g (^f p) (this matrix T is sometimes called the correlation matrix). The definitions 
for subsystem B are analogous. Since all p G C^) have = 0, we can encode any 

p G J^(C^^ C^) directly using only the remaining (M^ — 1)A^^ real values given by r"^ and T, as 
desired. 

Finally, we remark that in order to extend the result here to strong NP-hardness of EBP using 
the improved hardness results of Theorem [T| one must lower bound the distance /3' of T(<I>(p)) from 
the border of G2M,N, given the promise that p is P away from the border of Sm,n (to be clear, by 
distance here we mean the Frobenius norm of the difference of the Jamiolkowski representations 
of the maps). In particular, for inverse polynomial /?, we desire inverse polynomial in order to 
prove an inverse polynomial hardness gap for EBP. This may prove difficult, as T(<I>(p)) lives in a 
different space than p. We leave this as an open problem. 



Note that such a contained ball about I/MN does not exist when Qm^n is embedded into D(C^^ ® C'^) (as 
opposed to riC"' ® C'^)), since for any p e I?(C^'^ ® C^) with p ^ J^{C'^^ ® C^), and for any e > 0, we have that 
(1 - e)I/MN + ep is not in T{C'^' ® C^). 
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5 Concluding Comments 

We have shown that the problem of Weak Membership over the set of separable quantum states 
'Sm,n is strongly NP-hard, implying it is NP-hard even if one allows "moderate" error, i.e. /3 < 
1/ poly(M, A^). This shows that it is NP-hard to determine whether an arbitrary quantum state 
within an inverse polynomial distance from the separable set is entangled. It would be interesting 
to know whether WMEM^(5A/_Ar) remains NP-hard for /? a constant. Solving this problem would 
require a different approach than taken here, as, for example, even the first link of the reduction, 
CLIQUE <m RSDF, introduces an inverse dependence on the dimension. We also remark that the 
hardness result shown here is via a Turing reduction. It remains open, to the best of our knowledge, 
whether WM'EM./^(Sm,n) can be shown NP-hard under a many-one reduction. 

We next discussed two applications of our result: (1) We observed immediate lower bounds on 
the maximum Euclidean distance between a bound entangled state and Sm,n, and (2) we showed 
that the problem of determining whether a quantum channel is entanglement-breaking is NP-hard. 
Whether this latter problem is strongly NP-hard is left an an open question, for which we believe 
Theorem [6] should prove useful. 
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A Appendix 



Lemma 8. Given density operators pi,p2 G ^^(C 
P given by equations pi = + ■ a and p2 



*^®C^ 
I 

MN 



+ 



, with corresponding Block vectors a and 
\l3 ■ d (as per Eq. UV), respectively, we 



have 



\Pi - P2\\2 



1 

V2 



a 
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Proof. Via straightforward manipulation and the fact that Tr(crj(Tj) = 25ij, we have: 

1 



\Pl - P2\\2 



{oi - f3i)ai 



1=1 



Tr 



1=1 



M^N^-l 



V2 



12- 



□ 



Lemma 9. Combining Theorem]^ and LemmaVA gives ||c||2 G 0{\/N). 



Proof. By definition, we have: 



E 



-Tt(C<J,) 



(21) 



where m = M'^N'^ — 1. Recall now the definition of C from Eq. ([6]), where in our case, each 
Ai G R^^^ is all zeroes except for its upper left corner, which is set to submatrix Bi G H'^' 
from Theorem |2j Each Bi in turn is all zeroes, except for some index {k,l) (and hence {l,k), 
by symmetry), 1 < k < I < n, which is set to the (/c,/)th entry of the adjacency matrix Aq of 
graph G (see Theorem [2] and ensuing discussion) . We also require an explicit construction for the 
generators fij of SU{MN), given for example in [18], where {a^} 
that for 1 < p < g < MN and 1 < r < MN 
space n^'^^: 



iil^ ^ = {Upq,Vpq,Wr}, SUch 

1, and {xj}^^ an orthonormal basis for Hilbert 



Upq 
Vpg 

Wr 



ZXpXg -|- iXqXp 



rir + 1) 



rx,.+ix|;,^;L 



\k=l 



(22) 
(23) 

(24) 



Due to the symmetry of C and the fact that Tr(C) = 0, it is clear that only the generators of 



the form Upq will contribute to the sum in Eq. (21). Further, for each edge in G, Tr(C?7^ 



Jpq) 



for each Upq whose non-zero indices match those of the entries in C corresponding to that edge. 
Since each edge contributes four (symmetrically placed) entries of 1 to C, we hence have ||c||2 = 
^ (2e)22 = ^/2^, where e denotes the number of edges in G. Since e G 0{v?) {n the number of 
vertices in G), and N G 0(n^), we have ||c||2 G 0{^/N), as required. □ 
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